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Abstract 

We demonstrate the ability to measure the CP phase of a Higgs boson state by looking at the decays 
of heavy, Majorana neutrinos in the process pp — > H — > NN at the LHC. This is shown for signals 
with both opposite-sign (l^l^W^W^) and same-sign (l^W^W^) leptonic final states. These 
signals are investigated in the general framework of a sequential fourth generation of fermions with 
an additional right-handed neutrino. Such a scenario would naturally give rise to heavy Majorana 
neutrinos as well as significant enhancements to Higgs production via gluon fusion due to the 
contributions of the new, heavy quark states running in the fermion loops. Combined with the low 
background inherent to a same-sign lepton signal, this could be a useful way to investigate the CP 
nature of a Higgs boson state at the LHC. 
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I. INTRODUCTION 



While the Standard Model (SM) has been an extraordinarily successful theory of the 
strong and electroweak interactions, the precise origins of electroweak symmetry breaking 
are currently unknown. In the SM, this breaking is the result of the nonzero value for 
the vacuum expectation value (VEV) of a single scalar field, the Higgs boson, in a process 
known as the Higgs Mechanism. The SM Higgs boson has not yet been discovered, however. 
Furthermore, while the existence of a single Higgs boson - treated as a fundamental scalar 
particle - may be the simplest way to break electroweak symmetry, it is not the only way 
nor the most attractive; the weak-scale mass of the Higgs, and its stability under quantum 
corrections, are unexplained in the SM if it is to be understood as part of a larger theory with 
mass scale M 3> M\y, such as a Grand Unified Theory (GUT). In directly addressing these 
shortcomings of the SM, some theories such as weak-scale supersymmetry jl|, Little Higgs 
models Technicolor pf, Inert Higgs Doublet Models {4], and Twin-Higgs Models 5-7| 
feature an extended Higgs(-like) sector with multiple scalar bosons. Indeed, such extensions 
can be found in theories introduced to address other phenomenoloKical issues, including the 
origin of dark matter js-ll| and neutrino masses 
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A very interesting phenomenological consequence of multi-Higgs theories is the possible 
introduction of new sources of CP violation, which can arise from the complex mixing 
between the gauge and mass eig was, in fact, one of 

the primary motivations for first studying the Two Higgs Doublet Model (2HDM) 15|, Il6| . 
one of the simplest and most generic extensions to the SM Higgs sector. If a Higgs boson 
is discovered at the LHC or the Tevatron, determination of its CP nature will then be vital 
to more fully understanding electroweak symmetry breaking; detection of a CP-violating 
phase would, for instance, be an indication of an extended sector. 

A number of studies have discussed the detection of the CP phase of the Higgs boson 
using, for instance, angular observables sensitive to CP violation in the HZZ and HW + W~ 



couplings 17 
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27J as well as in the decays of the top quark in the process h — > tt — > bbW + W~ 



391 ] . The results in the latter case can apply equally well to other heavy fermions 



coupling to the Higgs boson and decaying via the W boson. Indeed, extensions of the SM 
with a fourth generation provide such heavy, fermionic candidates, a scenario which has 
received renewed interest (see, for e.g., Refs. 



40 



45J). 
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We consider here a fourth-generation model with a weak-scale Majorana neutrino. 
Through its mixing with the known fermions, it can decay to SM leptons via iV — > l^W^. 
The Majorana nature of these neutrinos then allows for signals with opposite-sign leptons 
and W bosons (H NN l+l'W + W~) as well as same-sign (H NN -»■ ^W^W^). 
The latter signal has long been understood to be the discovery mode for heavy Majorana 



neutrinos at colliders due to the low SM background (46|. We demonstrate here that both 
the opposite- and same-sign signals may be used in a manner akin to the top quark signals 
to determine the CP nature of the Higgs boson. 

In Sec. [Ill we describe a heavy, fourth-generation Majorana neutrino in a simple model 
and give its coupling to the SM Higgs boson and electroweak gauge bosons. We motivate 
the possible detection of this heavy neutrino in the gg — >• H — >■ NN mode at the LHC in 
Sec. IIIII We then discuss how, in the context of a more general model, the CP phase of 
a Higgs boson can, in principle, be determined via the angular correlations of the N decay 
products in the process pp — > H — > NN. 

II. FOURTH-GENERATION MAJORANA NEUTRINO 

We begin with a brief discussion of heavy Majorana neutrinos and their couplings to the 
Higgs boson in the context of a fourth generation 1 . We start with a simple model in which 
only the fourth-generation neutrino (N) develops mass, and therefore there is no mixing with 
the SM neutrinos. Consider the following Yukawa couplings and Majorana mass terms: 

Sf D —Y ri (f) T t Lr] u - X -M^ v + H.c, (1) 

where = (</> + , </>°) T is the usual SM Higgs doublet with hypercharge Y — 1/2, L — (x„, xi) T 
is the fourth-generation lepton doublet with hypercharge Y = —1/2, and ?7„ is a SM fermion 
singlet with hypercharge Y — 0. Here, all x an d rj fields are two-component, left-handed 
Weyl spinors and r = — ia 2 . After electroweak symmetry breaking, the neutral component 
of the Higgs doublet gets a VEV 



1 A more thorough discussion can be found in Ref. [471 ] and references therein. In this section we follow the 
conventions of Ref. 
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and this part of the Lagrangian becomes 



^D-Y n [v + -^Lj XuVu - \Mi) v T} v + H.c. (3) 



The mass terms can now be written as 



& 3 -\ (Xu, r)u) M n fa, Vuf + H.c. (4) 
with the Majorana mass matrix given by. 

(0 m D . 
(5) 
m D M 

Here we have defined the Dirac mass = Y n v. The mass matrix M n may be diagonalized 
by a unitary matrix U to give the mass eigenstates 



with the eigenvalues 
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When the Majorana mass term for the neutrino singlet, M, is much larger than we 
have Mi w m 2 D /M and M 2 w M. This is the usual see-saw mechanism invoked to explain 
the smallness of the masses of the known neutrinos when such a procedure is extended to 



the first three fermion generations (see Ref. 49] and references therein). For the following 



discussion we consider M to be large enough that N2 approximately decouples from our 
theory, leaving us with cos 6* ps 1 and one fourth-generation Majorana neutrino, TVi, that 
couples to the Higgs according to 

JSf D — j^HNxNi + H.c. (8) 
y/2v 

This can be written in the usual four-component notation with a Majorana spinor N = 
(N u Nl) T as 
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FIG. 1: Dominant pair-production subprocesses of a Majorana neutrino N through (a) the Drell- 
Yan process and gluon fusion to (b) Z and (b) Higgs bosons. 



JSf D -^HNP L N + H.c, (9) 

where we have relabeled the mass as Mn = M\. In this limit, N essentially acts as a gauge 
eigenstate when coupling to the W and Z bosons: 

2 D \-^-Z,N^P L N + i-Z= (W+N-fP L l 4 ~ W-l^P L N) , (10) 

a COS u\y y z 

where Z4 is the corresponding fourth-generation lepton whose mass is approximately ^> 
M^. The factor of i in the charged current, and relative minus sign between the two terms, 
is a remnant of the neutrino mixing: \ v = —icosONx —iNi. 



III. HEAVY NEUTRINO PRODUCTION: CP-CONSERVING CASE 

Let us now consider the pair production of such a heavy neutrino at the LHC in a model 
with a single, CP-conserving Higgs boson. The primary production modes are the Drell-Yan 
process and gluon fusion to the Z boson and Higgs boson via loops of heavy quarks. The 
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contributing diagrams are given in Figs. U(a)-(c). The corresponding partonic cross sections 
are given by 
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where v = y/2M w /g = (2V2G F y 1 / 2 w 174 GeV and 
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The normalized quark couplings to the Z boson are given by g\ = (T^ L — <5 9 sin 2 ^^) and 
g q R = (— Q q sin 2 6w), where T| L and Q q are the usual third component of weak isospin and 
electric charge, respectively, of the given quark. We note that both o z an d oh are a factor 
of 2 larger than the formulas in Ref. 51|, in which the pair production of heavy charged 
leptons were studied. This is due to the Majorana nature of the neutrinos we consider: a 
factor of 2 at the amplitude level from two possible Wick contractions of the neutrinos and 
a factor of 1/2 from identical particle phase space. 

We show the partonic cross sections in Eq. [TT1 convolved with CTEQ6L1 parton distri- 



bution functions 



50| for pp collisions at y/s = 14 TeV versus the neutrino mass for 



three different Higgs boson masses in Figs. |2] (a)-(c). We fix a = a(Mz) ~ 1/128 and use 
a running value of a s evaluated at /j,r = /jLf — 2Mjv, where //p is the factorization scale of 
our parton distribution functions and a s (Mz) = 0.118. For the range plotted in Figs. [2] (a)- 
(c), a s varies from ~ 0.1 to ~ 0.08. Here we have chosen degenerate fourth-generation 
quark masses of m ui = m d4 = mg 4 = 500 GeV. This choice, and our consideration of a 
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fourth-generation neutrino and Higgs boson masses of several hundred GeV, are consistent 
with experimental constraints as well as those from electroweak precision observables in the 



context of a fourth generation 42 1 . 



It is clear from Figs. [2] (a)- (c) that the Higgs boson contribution is dominant in all cases. 
This is largely due to the enhancement from the fourth-generation quarks 2 , a well-known 
enhancement of Higgs production from gluon fusion in the context of new heavy quarks 



51 



551 ] . In particular, in the range Mn « 100 — 200, when the process proceeds through 
a (nearly) on-shell Higgs boson, N pair production can have a cross section ranging from 
~ 100 fb to ~ 5000 fb. 

In Fig. [2] (d) we show, for direct comparison with the Majorana case of Fig. [2] (c), the 
pair production of a Dirac neutrino in a model with = 500 GeV. We note two distinct 
differences: the gluon fusion cross sections are a factor of 2 lower in the Dirac case, as 
discussed above, while the Drell-Yan process is actually enhanced. This latter effect is the 
result of a higher threshold suppression found in the Majorana case due to the axial-only 
coupling of the heavy neutrino to the Z boson. 

These results suggest that if a heavy neutrino related to a fourth generation is produced 
in pairs at the LHC, they may be directly associated with the Higgs boson (for M# ~ several 
hundred GeV) and they may be produced at detectable rates if they are not too massive. 
This link with the Higgs boson is particularly true in the Majorana case we study here, where 
the gluon fusion contributions are enhanced and the Drell-Yan contribution is su ppr essed. 



We note that this assessment agrees with those of previous studies, such as Ref. 



551 ] . who 



performed a thorough analysis of the rates of production and decay of heavy Dirac and 



Majorana neutrinos at the LHC. Similar scenarios have also been studied in Refs. 56l. 157]. 



2 When there is a large mass splitting (~ several hundred GeV) of these quarks, the gg —> Z —J- NN mode 
can be similarly enhanced and comparable to the Higgs contribution, a case studied in Ref. [51| for heavy 
charged leptons. However, corrections to electroweak precision parameters favor much smaller splittings, 



with ArriQ4 sa 50 GeV [42j . The production of heavy charged leptons from gluon fusion has since been 



revisited by Ref. [521 ] . who similarly find the Higgs contribution to be dominant. There can, however, be 



rather large contributions from an on-shell Z' decaying to NN, a case studied in Ref. [53 



7 




(c) (d) 

FIG. 2: LHC cross sections for pp — > NN, the pair production of a massive neutrino of mass 
M N , at yfs = 14 TeV with a Higgs boson mass of (a) 300 (b) 400 and (c) 500 GeV. The results 
for a Dirac neutrino with a Higgs mass of 500 GeV are given in (d). The short-dashed blue line 
gives the contribution from the Drell-Yan process qq — > Z — > NN , the long-dashed red line is for 
gg — > Z — > NN and the dotted-dashed magenta line gives the contribution from gg — > H — > AW. 
The thick black line gives the sum of the individual contributions. In all cases the fourth-generation 
quark masses are set to m„4 = = tuq^ = 500 GeV. 

IV. AZIMUTHAL DISTRIBUTIONS: CP- VIOLATING CASE 

Now that we have motivated the production rates of a heavy neutrino pair in a rather 
basic fourth-generation model, we consider two extensions. First, we allow small mixings 
with the first three generations 3 . This allows the heavier neutrino (s) to decay into SM 

3 See, for e.g., Ref. [5^ in which the four generations mix and there remains at least one neutrino of 
weak-scale mass along with three light neutrinos. 
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particles, which we parametrize with the neutral- and charged- current interactions 



& => l—^TT IKn ZpN-fPtVi + U UlN Z^PlN] (13) 

Zi COS ts'w 

+ [V lN W+N^PlI + v lN W-T>fP L N] , 

where v\ = u e , u^, u T , and I = e, fj,, r. Vin and U V[ n can, in principle, be different and we take 
each to be arbitrarily small but nonzero for at least one light lepton family. If, as described 
above, we take U and h to be heavier than N, the N — > l ± W T and N — > v\Z modes are the 
only relevant decays and the total decay width is very narrow. Along with very small values 
for Vin and U Ui n, this could lead to interesting effects at the LHC such as displaced vertices or 
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However, 



long-lived states which escape the detector and appear as large missing energy 
we assume that the widths are large enough [r^r > O(10~ 16 ) GeV] that the decays of the 
heavy neutrino are detectable inside an LHC detector. Also, independent of the absolute size 
of the mixing, we note that the maximum branching fractions to opposite-sign and same-sign 
final states are BR(l + l~W + W-) w 0.5 and [BR(l + l + W~W-) + BR(l-l~W+W + )] w 0.5, 
respectively, if there is only appreciable mixing to one generation and the mixing effects in 
the neutral current are suppressed. 

The second extension to the simple model presented in the previous sections, and the one 
of primary interest to this study, is the possibility of an extended Higgs sector with possible 
CP violation [ljj]. If, for example, we allow the Higgs boson state in Eq. [TUto be related 
to the mass eigenstates by a complex mixing element such that H = RuHi then Eq. [Ml 
becomes 

For a given Hi, the CP behavior depends on the precise nature of the mixing element Ru. 
For our purposes, we will parametrize the mixing in the Higgs sector, along with any residual 
mixing in the neutrino sector, by a complex parameter A to give the following coupling of 
the lightest Higgs to a heavy Majorana neutrino: 



if D -^HN(AP L + A*P R )N (15) 
D -^-HN(A R - iA n 5 )N, 

9 



where 



A=(A R + iAj) = \A\e 



(16) 



The phase of this coupling, a, is what we are ultimately interested in determining. Here, 
a = corresponds to a CP-even Higgs state and a = ±n/2 corresponds to a CP-odd 
Higgs state, where we restrict ourselves to —71/ 2 < a < it/ 2. Any nonzero value of a would 
indicate CP violation. In what follows, we seek to show that, in principle, the phase a may 
be determined by looking at the angular distributions 4 of the decay products of the heavy 
neutrinos in the process pp — >■ H — >■ iViV. In practice, the overall normalization of the cross 
section will be a mo del- dependent function of both the mixings in the Higgs and neutrino 
sectors as well as the precise nature of the Higgs' couplings to fermions. For example, in some 
2HDMs there are different Higgs bosons that give mass to the up- and down-type quarks. 
Here we take a mostly model-independent approach and present results that will depend 
only on Mn, Mh, and a and are independent of the normalization of a(pp — > NN). We 
also do not consider backgrounds, detector cuts and efficiencies, or the difficulties associated 
with reconstructing a boosted rest frame. 

A. Opposite-signs Case 

Given the above definitions, the squared matrix element for H — > NN — > l + W~l'~W + ', 
when summed over final-state spins, is 




(17) 



A] (2^1.(77,! + n 2 )u 2 .{n l + n 2 ) - (u 1 .u 2 )(n 1 + n 2 ) 2 ) 
+A\ (2u 1 .{n l - n 2 )u 2 .(rii - n 2 ) - (n 1 .n 2 )(n 1 - n 2 ) 2 ) 



+4A I A R e(u 1 ,u 2 ,n 1 ,n 2 ) , 



where the four-momenta are defined by 



4 These angular correlations were first studied in the CP-conserving case in Ref. 59 1. 
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m = l + + W~, (18) 

n 2 = + W + , 



Ul = ni -W (2-^ 

(™2-™2^ 



and Z ± , W ± represent the corresponding four-momenta of the leptons and W bosons. We 
have used the shorthand notation e(a, b, c, d) = e^ upa a^b v c p d a , where e is the antisymmetric 
Levi-Civita tensor with the signature e 0123 = +1. The factors of are related to the 
neutrino propagators and are given by 



(19) 



{n\ - MlY + (M N T N )i 
1 

712 | 



|2 



{nl-M 2 N y + (M N T N y 



The factor 6m is 1 if N is Majorana and is otherwise. We note here that this formula 
and the results of this section apply also for H — > tt — > bbW + W~ if we set 6m = and 
replace Mm — > M t , l~/l + — > b/b, and Vin — > V t b- In particular, Eq. [TT] matches the results 



of Ref. [l7j when rewriting their formula in terms of on-shell W bosons. 

We now define the z axis as the n\ direction in the Higgs rest frame and define the sets of 
angles (6*i,0i) and (6 2 ,4>2) of the W bosons (or equivalently the leptons) in the rest frames 
of the two heavy neutrinos using a coordinate system defined with respect to this axis 



\W X \ 

\W 2 \ 



sin 6*i cos 0i, sin 9\ sin <pi, cos ^i), (20) 
sin #2 cos 02, sin 9 2 sin <f> 2 , cos 9 2 ), 



where W\ = W and W 2 = W + . The relative azimuthal angle of the W bosons/leptons 
with respect to the n\ axis is $ = (02 — 0i)- This is the primary observable of interest and is 
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(a) 




(b) 

FIG. 3: A schematic of the decay H — > NN — > l + l~W + W~ in the rest frame of the Higgs boson, 
with a local coordinate system defined by z' = hi, the momentum direction of one of the heavy 
neutrinos. Note that z! is not, generally, aligned with the beam axis (typically denoted by z). The 
angle $ may be understood as the relative azimuthal angle of the decay planes of the two heavy 
neutrinos in the Higgs rest frame, as seen in (a). This is also demonstrated in (b), where W\p' 
and W2,T' denote the components of the VF-boson momenta in a plane normal to the z! = h\ axis 
and can be thought of as the transverse momenta of the W bosons in this rest frame (and are not 
to be confused with pt)- 



similar to that proposed in studies of CP correlations in h — > tt — > bbW + W (see, for e.g 



Ref. 



281 ]) and generic studies of CP- violating observables 601 ] . This geometry is illustrated 



in Figs.S(a)-(b). 

After making the narrow width approximation for and integrating over all angles except 
$, we obtain 
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— oc [ 16(M* + 2M 2 w f {PlA\ + A]) (21) 
-n 2 (M 2 N - 2M 2 V ) 2 (((3 2 N A 2 R - A 2 ) cos$ + 2(3 N A R A T sin®)} 

for the process H NN l + W^l'~W + . Here f3 N = - (AM 2 N /M 2 H ) is the velocity of 
the heavy neutrino in the Higgs rest frame, taken here to be on shell. As discussed above, 
the decay width of the neutrino should be extremely small, and therefore the use of the 
narrow width approximation here is very well justified 5 . 

The preceding equation can be written in a normalized form as 



1 dY 



1 

2^ 



7T 2 (M 



N 



2M 2 W ) 2 



{PnAr — Aj] 
lQ(M 2 N + 2M 2 w y W N A 2 R + A 2 ) 

Finally, we define a new angle x according to 



cos $ + 2 



PnArA! 



sin $ 



.(22) 



(/3 N A R ± iA T ) = (f3 2 N A 2 R + Aj)h ±ix 

that is related to a (the CP phase of the Higgs boson) by 

_ x / tanoT 



X = tan 



The $ distribution can now be given as 



V P. 



N 



(23) 



(24) 



1 dT 



i 

2^ 



M 2 

1 + 2^ 



cos($ - 2%) 



(25) 



The overall effect of the CP phase a in the coupling of the Higgs boson to NN is to 
introduce a phase shift in the cos($) dependence of the differential decay width. In the 
limit Mfj 3> 4M|r, we find x ~ * a an d the total phase shift is approximately 2a, providing 
a sensitive probe of the CP nature of the Higgs boson. In the other limit, Mjj « 4M^ and 
the NN are produced just above threshold. This leads to x ~ ±7r/2 for all but the smallest 



5 There are certain special cases, such as when « M\y, for which the narrow width approximation for 
./V may actually give large deviations from the true result even when (Tn/Mn) is small. We defer the 



reader to Rcf. 



6l| for further details. 
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values of a. This behavior is shown in Fig. H] (a). Measuring a in this regime would therefore 
require knowing (3n to very high accuracy. 

In general, the measurement of this phase will depend on the ability to detect the actual 
$ dependence of the distribution and, therefore, on the overall amplitude of the oscillation 
about the mean value. This is a function solely of the ratio My^/M^ and is plotted in 
Fig. IU(b). In the limit that Mn 3> Mw, the amplitude (relative to the mean value) reaches 
a maximum value of approximately (vr/4) 2 « 0.6. On the other hand, the amplitude vanishes 
for Mtv = \/2Mw- Therefore, the CP effects could be difficult to observe if N is not much 
heavier than the W boson. 

For completeness, we also give the result for H — > ti — > bbW + W~: 



1 dT 
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2^ 



where 



7T\ • 

4 ) 



9 Mk 



cos($ - 2 Xt 



(26) 



Xt = tan' 



tan a 



tan 



tana 



(27) 



Pt J \ y/l - 4M?/m 2 H 

This is the same as Eq. [25] with the replacement Mn — > M t . Indeed, for any heavy fermion 
of mass Mp coupling to the Higgs boson in a manner similar to Eq. [14] and decaying via 
F — > f'W , this formula will apply with the simple replacement M^r —> Mp independent of 
both BR{H -> FF) and BR(F -> f'W*). 



B. Same-sign Case 

While the opposite-sign signal has broad applicability to heavy fermions, we will now 
show a similar signal that is unique to Majorana fermions. The squared matrix element for 
H — y NN — > / ± /' ± H /=F W /=F , when summed over final-state spins, is given by 
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(71/2) (M w / M N ) 

(a) (b) 



FIG. 4: (a) The phase angle ^asa function of a for several different values of /3jv = y 1 — 4m^ / m? H . 
(b) The normalized, relative amplitude C of the $ dependence of the differential partial amplitude 
as a function of the ratio M W /M N . Here, C = (l - / (l + ' . 
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+(Wi o W 2 ). 
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Here we have ignored any interference with the crossed diagram, whose amplitude- squared 
contribution is represented here by the (W\ O W 2 ) term. This should be an excellent 
approximation, as the decay width of the Majorana neutrino should be extremely narrow. 
We can simplify the matrix element squared further if we only keep the terms that remain 
in the narrow width approximation: 

£Mf = -(2) 2 ( Vi^nw^ ) x |t ^ J2 x (29) 

|2«i.(ni + n 2 )u 2 .(n 1 + n 2 ) - (ui.u 2 )(rii + n 2) 2 
-4(^0^) (1 + ^)} 
+A^|2iti.(ni - n 2 )u 2 .(ni - n 2 ) - {ui.u 2 )(ni - n 2 ) 2 

+4(«i.ni)(« a .n a ) (l-^)} 

+4 A 7 A fl e(ui,M2, 711,712) 
+(Wi o W 2 ). 

There are some notable differences between Eq. [29] and the opposite-sign case in Eq. [171 
Apart from the {W\ -H- W 2 ) term, there are terms proportional to {ui.n\)(u 2 .n 2 ) that are 
absent from Eq. [TTJ In addition, the terms that Eqs. [29] and Eq. [T7] have in common differ 
by an overall minus sign. These differences ultimately lead to different behavior in the $ 
distribution. 

If we now label the leptons and W bosons, the four-momenta are given by 
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n 2 
m 

"2 



h + Wt, 
h + W 2 , 
n 1 -W 1 {2 



(30) 



7i 2 - Wa 2 



Ml 
(n 2 .n 2 ] 



The crossed term simply swaps the W boson labels in these definitions (e.g. rt\ — > l\ + W 2 ). 
In the discussion that follows, we can effectively ignore this term; it gives no contribution 
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to the overall decay width because of the compensating symmetry factor for the identical 
bosons, and any contribution to the angular distributions can be recast in terms of the first 
term by a relabeling of the particle states. 

If we now define $ and \ as before, we find that the azimuthal distribution for the W 
bosons/leptons is given by 

1 dT 1 

This differs from the result of Eq. [25] for the opposite-sign case only in the relative sign 
between the constant and cos($ — 2%) terms. It is therefore the same distribution shifted by 
an additional phase of ir. That the same-sign signal also shows a cos($ — 2x)-dependence 
in the differential partial width is of particular interest: this signal will have consider- 
ably less SM background than the opposite-sign case and will therefore be much easier 
to detect and reconstruct. In addition, when adding the H — > NN — > l + l + W~W~ and 
H — > NN — > l~l~W + W + contributions together, the overall rate should be the same as 
in the opposite-sign case. The same-sign signal would therefore be much better suited for 
making a determination of the factor 2\ and, hence, a and the CP nature of the Higgs 
boson. 



7T 



M 2 

1U N 



cos($ - 2%) 



(31) 



V. COMPARISON TO MADGRAPH 

We perform a simple check of our results by implementing a model with a heavy Ma- 
jorana neutrino (with cou plin gs to the Higgs and W boson as described above) in Mad- 



Graph/MadEvent v4.4.57 |62j. To estimate its decay width, we assume the heavy neutrino 
decays entirely through W bosons with Vw ~ 10~ 3 , such that ~ 2|V/7v| 2 r t ~ 10~ 6 GeV 
(where ^ « 1.5 GeV is the width of the top quark). We look at the decays H — > NN — > 
l + l~W + W~ and H ->■ NN l ± l ± W T W T and in each case generate 50000 events in the 
rest frame of the Higgs boson. In Figs. [5] (a)-(d) we present the partial decay width as a 
function of $ for the opposite-sign signal using Mh = 500 GeV and = 175 GeV. This 
is plotted for a CP-even Higgs boson (a = 0), a CP-odd Higgs boson (a = 7r/2), and 
mixed-CP state (a = vr/4). In all cases, the MadGraph results show the expected behavior. 
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The same results are given for the same-sign signal in Figs. [6] (a)-(d). The analytic curves 
match the MadGraph results quite well and demonstrate the overall phase difference of it 
with the opposite-sign results, as described above. 

In both cases, it is important to properly pair the leptons and W bosons to reconstruct two 
on-shell neutrinos; in the opposite-sign case, the charge of the W bosons may be ambiguous 
if they both decay hadronically and in the same-sign case there is a fundamental ambiguity 
even if the charges are known. We followed a very simple procedure here: as we know each 
event is derived from the decay of two equal mass intermediate states with a very narrow 
width, we take the combination that minimizes + Wi) 2 — (I2 + Wj) 2 |. This works rather 
well for our MadGraph-level results; however, due to smearing effects and the inherent 
difficulties that will arise in reconstructing the Higgs rest frame, more elaborate techniques 
will likely be necessary in practice. 

In Figs. [7] (a)-(d) we show, for comparison, results derived from the opposite-sign signal 
for a heavy neutrino with mass just below production threshold, = 240 GeV. This 
demonstrates both the enhancement to the relative amplitude due to the small value of 
Mw/Mn as well as a phase-shift 2\ ~ 7r in the case of a 7^ with near-threshold production. 

VI. CONCLUSIONS 

We have demonstrated that the CP phase of a Higgs boson state in an extended Higgs 
sector may, in principle, be determined at the LHC by looking at angular correlations in the 
signals H -> NN l + l~W + W~ and H -> NN -> l^W^W*. This is achieved by looking 
at the differential partial decay width as a function of $, the relative azimuthal angle of the 
leptons / W bosons about the axis defined by the neutrino momenta in the Higgs rest frame. 
The CP phase of the Higgs boson, a, introduces a phase-shift 2\ = 2tan _1 (tan(a;)//37v) i n 
the cos<3> dependence of the decay width, where (3n is the velocity of the neutrinos in the 
Higgs rest frame. Measurements of this phase shift will allow a direct determination of the 
CP nature of the Higgs boson with this signal. 

We have also verified that in the context of a fourth generation with heavy quarks, the 
process pp — > H — > NN can have a cross section of ~ 100 fb - 5000 fb at the LHC if 
the Higgs boson is heavy enough for on-shell decays to heavy neutrinos with mass > 
100 GeV. This process, including the subsequent N decays, may therefore be observable 
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9 1 1 1 1 ■ 1 1 1 1 ■ 1 o 2 1 1 ' 1 1 1 ' ' 1 1 1 
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(*/7l) (<I>/71) 

(c) (d) 

FIG. 5: Angular distribution of <5 for opposite-sign leptons/VF bosons using Mh = 500 GeV, 
M]\r = 175 GeV, and a CP phase of (a) a = 0, (b) a = f , and (c)q = \. The solid curves give 
MadGraph results and the dashed curves are our analytic predictions. The theory results for the 
three choices of phase are plotted together in (d) . 

above background, particularly in the case of the same-sign signal, which may allow for the 
azimuthal distributions to be measured at the LHC. 
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FIG. 6: Angular distribution of $ for same-sign leptons/iy bosons using Mh = 500 GeV, Mn = 
175 GeV, and a CP phase of (a) a = 0, (b) a = ^, and (c)a = f . The solid curves give MadGraph 
results and the dashed curves are our analytic predictions. The theory results for the three choices 
of phase are plotted together in (d). 
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